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ABSTRACT 

This paper presents an analytical procedure to simulate vibrations in 
gear transmission systems. This procedure couples the dynamics of the rotor- 
bearing gear system with the vibration in the gearbox structure. The modal 
synthesis method is used in solving the overall dynamics of the system, and a 
variable time-stepping integration scheme is used in evaluating the global 
transient vibration of the system. Locally each gear stage is modelled as a 
multimass rotor-bearing system using a discrete model. The modal 
characteristics are calculated using the matrix-transfer technique. The 
gearbox structure is represented by a finite element model, and modal 
parameters are solved by suing NASTRAN. The rotor-gear stages are coupled 
through nonlinear compliance in the gear mesh while the gearbox structure is 
coupled through the bearing supports of the rotor system. Transient and 
steady state vibrations of the coupled system are examined in both time and 
frequency domains. A typical three-geared system is used as an example for 
demonstration of the developed procedure. 
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Torsional damping matrix 

Bearing excitation forces 

External excitation forces 

External excitation moment 

Gear mesh force in x- and y-directions 

Gear mesh torque 

Gyroscopic-angular rotation matrix 
Gyroscopic-angular acceleration matrix 
Identity matrix 

Rotational mass moment of inertia matrix 
Average stiffness matrix 
Gearbox stiffness matrix 

Compensation matrices in x- and y-direction 
Gear mesh stiffness between i th and k th rotor 
Shaft stiffness matrix 

Bearing direct and cross-coupling stiffness matrix 

Torsional stiffness matrix 

Mass-inertia matrix of rotor 

Mass-inertia matrix of gearbox 

Radius of gear in the i th rotor 

Gear generated torque 

Generalized motion in x- and y-directions 
Gearbox motion in x-, y- and z-directions 
Gearbox motion at bearing supports 
Gear displacements in x- and y-directions of 
the i th rotor 

Gear forces in x- and y-directions 


2 


V V 

a ki 
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Motion of rotor at bearing support 

Angle of tooth mesh between k th and i th rotor 

Rotational displacement of the i th gear 

Lateral and torsional eigenvalue diagonal matrices 

Lateral and torsional orthonormal eigenvector 

matrices of the k th rotor 

j th orthonormal mode of k th stage at 1 th node 


INTRODUCTION 

Recently there has been an increase in the use of gear transmissions in 
both defense and commercial applications. The ever increasing speed and 
torque requirements of newer transmission systems , often result in excessive 
noise and vibration at both the gear stages and the gearbox structure. Today 
there is a wealth of literature concerning noise and vibration reduction 
through gear tooth modification and design. However, the study of the 
dynamics and acoustics of the overall gear transmission system is somewhat 
limited. 

August (1986) studied gear system vibrations for a planetary gear system. 
Boyd and Pike (1987) and Choy (1988b) used the work done by Cornell (1981) to 
study the dynamics of various gear transmission systems. Mitchell (1985) and 
David (1987, 1988) simulated the dynamics of multistage gear systems using the 
matrix transfer method. Choy et al. (1989, 1990) calculated the dynamics of 
the multistage gear systems with effects of base motion using the modal 
method. Ozguven and Houser (1988) and Kahraman et al. (1990) used a finite 
element model to predict the dynamics of multistage gear systems. In term of 
gearbox vibration analysis, some work has been reported by Lim (1990) using 
finite element analysis. Very little work, however, has been cited in the 
literature concerning gearbox coupled vibration in gear transmission systems. 
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The work presented in this paper is the development and application of a 
combined approach of using the modal synthesis and finite element method in 
analyzing the dynamics of multistage gear systems coupled with the gearbox 
structure. Modal equations of motion are developed for each rotor-bearing- 
gear stage, using the matrix transfer method, to evaluate the modal 
parameters. The modal characteristics of the gearbox structure are evaluated 
using a finite element model on NASTRAN. The modal equations for each rotor 
stage and the gearbox structure are coupled through bearing supports and gear 
meshings. The modal equations are solved simultaneously with the appropriate 
initial conditions. The modal accelerations are integrated using a variable 
time-stepping integration scheme to obtain the transient vibration of the 
system. A typical three stage gear system is used as an example for this 
analysis. Results are presented in both time and frequency domain and in both 
modal and generalized coordinates to facilitate a complete representation of 
the dynamic characteristics of the system. 

Development of Equations of Motion 

The equations of motion for a single stage multi-mass rotor-bearing-gear 
system with the coupling effects of gear-box vibrations and the rotor inertia- 
gyroscopic effects can be written in matrix form for the i th stage (Choy, 

1987; 1989) for the X-Z plane as: 

[M] 1 {X} i + lG v ) i {iy 1 + [c xx ] 1 {i -*„>! + - Vi 

+ [ g aM y >i + t R xx *k 8 ] i{ x }i - r K xxM x bs>i 

+ ( K xyM Y - Y bs>i " <*«<*>>! + < F G x(t)i> U> 
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and in the Y-Z plane as; 

[M]i{?>i - [G v ] i {i} i + [C yx ] i{ £ - X b } i + [G yy ] i {^ “ Y b } i 

“ C G A ]l< X >i + t K yy +K sli< Y >i - [ K yyh< y bs>i 

+ [K yx ] i{ Y - Y bs}l - {F y (t)}. + {F Gy (t) }i (2) 

Here F and F are force excitations from the effects of mass imbalance 
X y 

and shaft residual bow in both X- and Y-directions . F Qx and F Qy are the X 
and Y gear mesh forces induced from the gear teeth interaction with other 
coupled gear stages. The bearing forces are evaluated through the relative 
motion between the rotor {X}, {Y} and the gearbox {X b >, {Y fa } at the bearing 
locations (Choy, 1987). The mass-inertia and gyroscopic effects are 
incorporated in the mass matrix [M] and the gyroscopic matrices [G y ] and [GJ. 
The coupled torsional equations of motion for the single rotor-bearing-gear 
system can be written as: 

[JJitB > 4 + [C T ) i {8 > t + [K, = {F T (t)} i + {F ot (t)> 4 < 3 ) 

In Eq. 3, {F (t)} represents the externally applied torque and {F Qt ( t ) } 
represents the gear mesh induced moment. Note that Eqs. (1) to (3) repeat for 
each single gear/rotor stage. The gear mesh forces couple the force equations 
of each stage to each other as well as the torsional equations to the lateral 
equations (Choy, 1989; Cornell, 1981; and David, 1987; 1988). The coupling 
relationships between the torsional and the lateral vibrations and the 
dynamics of each individual gear/rotor stage are derived in the next section. 
In addition, there are equations of motion for the gearbox which couple the 
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various rotor stages through the bearing supports. The gearbox equations can 
be written as: 

X-equation 

[M b H* b > + [C bx ]{* b > + [C xx ]{i b -* 8 > + [C xy ]{* b -t„> * [K bx ]{x b > 

- [K xx ){X b -X 8 } + [K xy ]{Y b - Y s } » 0 (4) 

Y-equation 

(M b ]{? b > + [c by ]{* b > + [C yx ]{i b -i a > + [C yy ]{i b -i 8 > + [K by ]{Y b > 

- f K yyHV b -Y s > + [K yx ]{X b -X 8 > =0 (5) 

and Z-equation 

[M b ]{2 b } + [C bz ] {t b y * [K bz J{Z b } -F bz (t) (6) 

where F (t) is an excitation function due to external forces in the axial 

bz ' 

direction . Since the bearing is assumed to be uncoupled in the Z-direction, 
Eq. (6) can be solved independently without considering shaft motion. 

Coupling of Gear Meshes 

The torsional and lateral vibrations of a single individual rotor and the 
dynamic relationships between each gear stage are coupled through the 
nonlinear interactions in the gear mesh. Gear mesh forces and moments are 
evaluated as functions of relative motion and rotational between two meshing 
gears and the corresponding gear mesh stiffnesses. These gear mesh 
stiffnesses vary in a repeating nonlinear pattern with each tooth pass 
engagement period (August , 1986 and Cornell , 1981) and can be represented by a 
high order polynomial (Cornell, 1981 and Boyd, 1987). A sixth order 
polynomial curve is used in this study to simulate the stiffness changes for 
contacting gear pairs (zero stiffness are input for noncontacting pairs). The 
repeatability of such nonlinear mesh stiffnesses can also act as a source of 
steady state type of excitation to the gear system. With the coordinate 
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system as shown in Fig. 1, the following gear mesh coupling equations can be 
established by equating forces and moments (Choy, 1989). For the k th stage 
gear of the system, summing forces in the X-direction results in: 

n 

F Gxk ’ E K tKil- R ci 8 ci- R cK 9 ck + ( X ci " X C k> COS fl ki 
1=1,1 *k 

+ < Y ci - Y ck) sin a kil COS a ki (7 

Summing forces in the Y-direction results in: 

n 

F Gyk - £ K tki[- R ci 0 oi- R ck 0 ck + < X ci - X C k> COS «ki 

1=1,1 

+ (Y cl - y ck )sin a kl ]sin a kl (8 

Summing moments in the Z-direction results in: 

n 

F Gyk " E R ck {K tkil { -R ci 0 ci -R C k 0 ck) + ( X ci - X ck) COS «ki 
1=1,1 7*k 

+ (Y cl - Y ck )sin a ki ]> (9 

where n is the number of stages in the system. 

Modal Analysis 

To reduce the computational effort, the number of degrees-of-f reedom of 
the system are reduced through modal transformation. Orthonormal modes for 
each individual rotor-bearing stage are obtained by solving the uncoupled 
system homogeneous characterist ic equations. Using the modal expansion 
approach (Choy, 1987; 1988a; and 1989), the motion of the system can be 
expressed as: 


m m 
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E A bi< 0 bxi> 
1-1 
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where m is the number of modes used to define each motion. The 
orthogonality conditions of the modes can be expressed as: 

[$] t [K] [*] = [A 2 ] 

where 

K - < K xx + K yy )/ 2 + K s 

[$ b ] T [K b ] [*„] = [A*] 

[* t ] T [K T ] [* t ] - [Aj (13) 

[*] T [M][$] = ($ b ] T [M b ] [$ b J - [I] (14) 

Using the modal expansion and the orthogonality conditions, with the 
bearing forces due to the base motion expressed in the right hand side of the 
equation, the modal equations of motion for the rotor bearing system 
(Choy, 1989) can be written as: 

X-Z equation 

{A> + [$] t [g v ] [$]{6> + [$] T [C xx ] [$]{A> + ($] T tc xy H$]{6> 

+ [$] T [G a ][$]{B} + [A 2 ] {A} + [#] T [K dx ][$]{A} 

+ {$] T [K xy ][$]{B} - ($] T {F x (t) +F Gx (t) + F Bx (T)} (15) 

where 

t K dx] - t K xx) + r K a ] - < 16 > 

and. 


(Hi 


( 12 ) 


F Bx (t) - [C xx ]{Z b ) * [C xy ]{^ b > + (K xx ]{X b > + [K xy ){Y b > (17) 
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which can also be expanded into modal parameters as 


F Bx( fc ) = l c xxH* bx ]{A b } + [C xy ][$ by ]{6 b } 

+ l K xxH* bx HA b > + [K xy ][$ by ]{B b > 

Y-Z equation 

{6} - [$] t [g v ]($]{A> + [$] t [c yx ][$]{A> + [♦] T [C yy ][f]{6> 

-[♦] T [G a ] 1*]{A} + (A W J {B} + [*] T [K dy ]($]{B} 

♦ [^] T [Kyx][^]{A> = [^] T {F y (t) +F 0y (t) +F By (t)} 

where 


[K^] = [K^] + [ Kg ] - [K] 


F By (t) - [C yX ]{i b } + [C yy ]{^} ♦ [K yx HX b } + [Kyy){Y b } 


which can also be expanded into modal parameters as 


F B y<t) - [C vx ][$ bx ]{A} + [C vv ][$ bv ]{6} + [K vx ][^ bx ]{A> + [K yy ][$ by ]{B} 


“'yx J i 


“'yy J 1 by 


yx J LTbx J 


and the 8-equation can be expressed as 

{A t > + [$ t ] T [C t ][$ t HA t > + [Aj{A t > = [* t ] T {F t (t) +F ot (t)> 


The gear mesh forces and moments can also be expressed in the modal 
for the k th stage with the gear location at the 1 th node, as: 


] k {F Gx > - £ 0k jl 

3-1 


E K t*i[ ~ R ci 0 ci ~ R c k d c k + < X ci - X ck> COS a 
l-i, iA 


kl 


+ < Y ci - Y ck> Sin “kil COS a ki 
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( 23 ) 
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I*M F Gy> * £ 0k 3 L 


J-l 


E " R ci 0 ci — R ck ^ck * < X cl " X ck> COS «. 


ki 


i*i # iA 


+ < Y ci - Y ck> Sin “kJ Sin a 


( 25 ) 


ki 


t*M F ot> * £ a ^ 1 


j-i 


£ R ck( K tkif < -*cAi - R ck 0 ck) + < X ci - X ck> 

i=l,iA 


+ cos a kl + (Y ci - y ck ) sin a ki ] ) 


(26) 


where k is the stage number, j is the mode number, and 1 is the station 
location of the gear mesh. 

A set of modal equations of motion can also be written for the gearbox 

as: 

X-equation 


l*bx] T [ M bH* b xHA b > + [* bx ] T tc bx ][$ bx ]{A b > + [* bx ] T [K bx ][$ b J{A b } 

+ t $ bxl T <t K xxHX b -x 9 > + [C xx ]{x b -* 8 > 

+ fK xy 3{Y b - Y s > + [C xy H* b - * a >> - 0 (27) 


For proportional damping in the gearbox model, the equation can further 
be reduced to: 


[i}{A b ) + [C x ]{A b > + [A b ] {A b > + [^ bx ] T {F Bx (t )} - [* bx ] T {[c x ][*){A> 

+ [C xy )($){6> + (K XX ][$]{A> + [K xy H$]{B}> =0 (28) 
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Similarly/ the Y-equation can be written as 


[I]{fi b } + [C by ]{6 b > + [Aj{B b > + [$ by ] T {F By (t)} - [$ by ] T {[c yx n$]{A> 

+ [C yy ][$){6} + [K yx ][$]{A} + [K yy ][$]{B>} - 0 (29) 

and the Z-equation as 

[ I] <0 b > + l C bzHV + [A]{D b > = [* bz ] T {F bz (t) } (30) 

Solution Procedure 

In order to obtain a solution for the overall dynamics of the gear 
transmission system/ a three phase solution procedure is used, namely, 

(1) (the evaluation of modal characteristics of the rotor-bearing systems and 
the gearbox structure, (2) transient vibration solution of the overall 
dynamics in modal coordinates, and (3) the transformation of force and 
vibration data from modal into generalized coordinates. Discussions of the 
three solution phases are presented in following paragraphs. 

The natural frequencies and mode shapes required in the transformation of 
the equations of motion for the rotor-bearing system (Eqs. (1) to (3)) into 
the modal equations (Eqs. (15), (19), and (23)) are evaluated using the matrix 
transfer method (Choy, 1989 and 1990) on a discretized lumped mass model. 

Using an average stiffness value for the bearing supports, the undamped modes 
are calculated using a marching search technique for the assigned frequency 
range. The natural frequencies and mode shapes for the gearbox structure are 
obtained by the finite element analysis using NASTRAN. The modal data from 
the gearbox analysis are used to transform the gearbox equations of motion 
(Eqs. (4) to (6)) into their modal coordinates (Eqs. (28) to (30)). 

The second phase of the solution procedure involves the solution of the 
coupled modal equations of motion between the rotor stages and the gearbox 
structure. The coupling effects of the gear mesh and the bearing supports are 


11 


also expressed in modal coordinates such that the global equations are solved 
simultaneously in modal form* Using a set of initial conditions for both 
displacement and velocity of the global system, calculated from the steady 
state conditions at the rotor-bearing systems and zero vibration at the 
gearbox, the modal accelerations A, B, A fa , B^, and of the system can 

be evaluated (Eqs* (15), (19), (23), and (28) to (30)). A variable time 

stepping integration scheme (the Newmark-Beta Method) is used to integrate the 
modal acceleration to evaluate the modal velocities and displacements at the 
next time step. A regular time interval of 200 points per shaft revolution is 
used in this study with a refined region of smaller steps at the gear mesh 
transition period for single and multiple teeth contact. 

The modal acceleration, velocity, and displacement calculated from the 
transient integration scheme are transformed back into the generalized 
coordinates (Eq. (10)). The nonlinear bearing forces and gear mesh forces can 
be evaluated from the velocity and displacement differentials between the 
rotor stages and the gearbox structure. Results from this solution procedure 
are demonstrated in the following section using a prototypical gear 
transmission system. 

Discussion of Results 

To demonstrate the application of the discussed analytical approach, a 
three-gear transmission system given in Fig. 2 is used as an example. Note 
that the gearbox is assumed to have uniform thickness throughout the enclosed 
walls and is fixed to the ground at the four lower corners. Figure 3 shows 
the arrangement and orientation of the gear stages inside the gearbox 
structure. While all three gears have an identical 36-tooth gear and a mesh 
contact ratio of 1.6, the driver (stage 1) is longer and larger in diameter 
than the two other stages, which are identical in geometry. Stage 1 is 
supported by two bearings located at the end plates while gear stages 2 and 3 
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are supported by bearings located at both end plates and at the middle of the 
gearbox, as shown in Fig. 3. An operating speed of 3000 rpm is used for all 
three gears. 

In order to calculate the global dynamics of the gearbox system, the 
modal characteristics for both the gearbox structure and the rotor-bearing- 
gear stages are evaluated. The finite element approach (NASTRAN) is used to 
model the gearbox structure using a combination of plate elements. Table I 
presents the results of the first nine natural frequencies of the gearbox 
structure with their corresponding three-dimensional mode shapes given in 
Figs. 4 to 6. The dynamics of each individual rotor-bearing-gear stage are 
modelled using the matrix-transfer method. Using an averaged bearing 
stiffness, as discussed in the previous sections, the lateral natural 
frequencies and mode shapes of each rotor stages are evaluated. A similar 
procedure is repeated to evaluate the torsional vibrations of the rotor 
stages. Table II presents the results of both torsional and lateral natural 
frequencies for all three rotor stages. Some of the lateral and torsional 
mode shapes for the rotor stages are given in Figs. 7(b), 8(b) and 9(b). 

Using mass imbalance in all three rotor stages as excitation input 
and a nonlinear gear mesh compliance between the gear stages, the global 
transient dynamics of the system are evaluated using zero initial conditions. 
Figures 10 to 13 present the gearbox vibrations in terms of modal excitations 
of the first 8 natural frequencies. These modal excitations (Choy, 1987; 
1988a; and 1989) represent the excitability of the particular mode and are 
expressed in the frequency domain using an FFT procedure to transform the 
modal time variables into the frequency domain. Note that the major component 
excited in each mode is at its own natural frequency. The highest component 
of excitation is in the x-direction. A moderate excitation is seen in the 
y-direction, while a very small magnitude of excitation exists in the 
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z -direction. A closer examination shows that the two highest x components 
occur at the second and fifth mode (462 and 575 Hz) while the highest in the 
y-direction occurs at the third mode (509 Hz). Figure 14 presents the total 
vibration at the upper corner of the gearbox (node 82) in both time and 
frequency domains. Note that three dominant components excited in the 
x-direction are located at frequencies of 462 , 509, and 575 Hz with the 
highest two located at 462 and 575 Hz. The Y-direction has only one major 
component a 509 Hz while no significant vibration is detected in the 
z-direction. This further confirms the use of modal excitation parameters 
in representing the dynamic behaviors of the global system. 

Figure 15 shows the orbiting motion of the gear stages during the initial 
transient period. Note that the vibration of the system eventually settles 
into a steady state motion. The smaller first stage orbit is due to the large 
stiffness in the rotor-bearing system. The gear mesh forces between gears in 
stages 1-2 and 1-3 in both time and frequency domains are given in Fig. 16. 
Three major frequency components exist in the gear mesh forces at 0, 50, and 
1800 Hz. They represent the static gear load at 0 Hz, the vibratory frequency 
at rotational speed of 50 Hz, and the tooth pass frequency at 1800 Hz. The 
existence of sidebands at 1800 Hz is mainly because of the shifts in the 
toothpassing frequency due to the vibration of the rotor-gear system. A 
further amplification of such frequency components can be seen in Figs. 7(a), 
8(a), and 9(a). The modal excitation of the first lateral mode in Fig. 7(a) 
has its largest component at 50 Hz due to mass imbalance at rotational speed. 
The modal components of the first torsional mode in Fig. 8(a) show that the 
gear load at zero frequency and the tooth pass frequency vibration are excited 
in stages 2 and 3. Figure 9(a) shows that the zero gear load frequency is 
excited in stage 1 at the second mode. This again further confirms the use 


14 


of the modal synthesis approach in analyzing dynamic behaviors in gear 
transmission systems, 

SUMMARY 

This paper presents a vibration analysis of a gear transmission system in 
which the dynamics of the gearbox structure is coupled with the vibration of 
the rotor-bearing-gear stages. The analysis combines the modal 
characteristics of the gear stages developed through the matrix transfer 
method with the modal parameters of the gearbox evaluated by a finite element 
model. The major content of this work can be summarized as follows: 

1. A comprehensive procedure is developed to combine the dynamics of the 
rotor-gear system with vibrations of the gearbox structure to determine the 
global system response. 

2. The modal method is used to transform the equations of motion into 
modal coordinates before the synthesis of the global systems of equations of 
motion to reduce the degrees of freedom of the system. 

3. The use of modal excitation functions in the frequency domain provide 
good insights of the dynamic behavior of the gear transmission system. Such 
knowledge is crucial for designing transmissions with improved performance and 
durability . 

4. The sensitivity of the gearbox vibration, due to rotor mass imbalance 
and gear mesh nonlinearity, can be evaluated using the developed methodology. 

5. The initial transient dynamic analysis approach developed in this 
study can also be applied to simulate conditions in which a sudden excitation 
is applied to the gear transmission system. 
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TABLE I. - GEARBOX 
NATURAL 
FREQUENCIES 


Mode 

Hz 

1 

412 

2 

462 

3 

509 

4 

540 

5 

575 

6 

697 

7 

852 

8 

893 

9 

911 


TABLE II. - ROTOR-GEAR STAGE NATURAL FREQUENCIES 



Lateral mode, 
Hz 

Torsional mode, 
Hz 

Axial mode, 
Hz 

Stage 1 

131 

661 

3 931 


297 

2086 

5 208 


1047 

4197 

11 724 

Stage 2 

123 

186 

1 918 


187 

866 

4 860 


390 

1686 

8 186 

Stage 3 

123 

186 

1 918 


187 

866 

4 860 


390 

1686 

8 186 
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Figure 5. — Gearbox 3-D mode shapes for modes 4, 5, and 6 
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Figure 8. — Rotor first torsional modal excitation and mode shapes for stages 1, 2, and 3. 
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Figure 10. — Gearbox modal excitation for modes 1 and 2. 
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Figure 1 1 . — Gearbox modal 
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Figure 12 —Gearbox modal excitation for modes 5 and 6. 
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